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Abstract 

The position of an interface (domain wall) in a medium with random pinning 
defects is not determined unambiguously by a current value of the driving force even 
in average. Based on general theory of the interface motion in a random medium 
we study this hysteresis, different possible shapes of domain walls and dynamical 
phase transitions between them. Several principal characteristics of the hysteresis, 
including the coercive force and the curves of dynamical phase transitions obey 
scaling laws and display a critical behavior in a vicinity of the mobility threshold. At 
finite temperature the threshold is smeared and a new range of thermally activated 
hysteresis appears. At a finite frequency of the driving force there exists a range 
of the non-adiabatic regime, in which not only the position, but also the average 
velocity of the domain wall displays hysteresis. 

Key words: 
PACS: 



1 Introcuction 



This is our tribute to the memory of Per Bak, great scientist and an excep- 
tional personality. His main contribution to science is the discovery of a new, 
very broad class of nonlinear deterministic systems displaying stable chaotic 
critical behavior, which he called Self- Organized Criticality (SOC). It is diffi- 
cult to overestimate the significance of this discovery, which has opened new 
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ways for explanation and unified description of very remote phenomena, such 
as earthquakes, plastic deformations, hysteresis in magnets, dynamics of bio- 
sphere and financial catastrophes. The introduced by Per notion of avalanches 
is now central for dynamics of such systems and is most probably a main 
source of the 1// noise. Theory of these systems has deserved a wide popu- 
larity and Per is one of the most cited physicists. 

However, his moral influence in the scientific community was not less impor- 
tant than his scientific contribution. He always was surrounded by people seek- 
ing new ways in different fields: geologists, biologists, physicians, economists, 
historians and, certainly physicists. Everybody felt himself free in this com- 
munity to criticize everything and to propose new ideas, but Per's authority 
and respect to him was extremely high. It was based on his uncompromising 
pursuing of truth and hate to any falsehood in science as well as in social 
and personal fife. Once he was shown by a big TV company at JFK airport 
expressing his indignation by the violation of passenger rights. He was a brave 
and noble man. He could be very harsh in discussions, a feature, which not 
everybody could tolerate. However, Per was a reliable, friendly and responsive 
person ready to help if there was a need in his help, and for this purpose he 
could spend his time, money and infiuence. We feel that his death has left a 
vacancy in our community, which will not be filled. 

In this review we consider a simple, but rather wide-spread phenomenon: the 
hysteresis mediated by a motion of an interface or domain wall (DW) driven 
by an external alternating force in a medium with random pinning centers. 
The interplay between the driving force and elastic and pinning forces devel- 
ops in time. If the driving force is constant, it can establish an average velocity 
of the DW after a sufficiently long time interval. In the case of the alternating 
force this time may be longer or shorter than period of oscillations resulting 
in different regimes. The infiuence of a heat reservoir gives an additional di- 
mension to this phenomenon. The interface is a typical system displaying the 
S0C[1]: in the absence of the external force it exhibits a fractional structure. 
A small force results in avalanches on spatial and temporal scales depending 
on the force. The stronger is the force and the higher is the frequency the 
weaker is the infiuence of the avalanches, but they are crucial for the range of 
small frequencies and amplitudes important in the experiment. 
This problem is a part of a very old hysteresis problem [2] occurring in many 
different dynamical processes: chemical reactions, magnetization reversal, crys- 
tal growth, absorption and desorption etc., in which many DW are involved. 
In this case the interaction between walls must be taken in account. It is be- 
yond the scope of our consideration, but some rough estimates can be done 
on the basis of one-DW theory. On the other hand, the recent experimental 
development with its tendency to sub-micron spatial scale generated many 
small systems in which the motion of one or few DW is absolutely realistic. 
Moreover, the motion of an individual domain wall was experimentally dis- 
covered and studied by several groups of experimenters [3]. Therefore, theory 
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of such a motion is necessary for understanding of the observed phenomena 
and prediction of new ones. 

Though the observation of the hysteresis in the ferromagnets accounts more 
than 100 years history, the phenomenon could not be satisfactory explained 
until the end of 20-th century, when a theory of DW and its motion in a 
random medium was developed. Therefore we start our review with an intro- 
ductory section briefly describing this theory. In the next section we consider 
a DW moving adiabatically at zero temperature. Adiabatic motion means 
that the frequency of the driving force is small enough and the momentarily 
velocity is equal to its stationary value at a constant force equal to the mo- 
mentarily value of the ac force. The third section is dedicated to the influence 
of flnite temperature, but still in adiabatic regime. In the fourth section we 
consider non-adiabatic regime predict a new phenomenon: the hysteresis of 
velocity (magnetization rate in the case of a ferromagnet). In the last section 
we review relevant experiments. 



2 Domain wall in a random medium 

2.1 Zero temperature 

For deflniteness we consider a DW in an impure ferromagnet with uniaxial 
anisotropy at zero temperature. As it was shown in [4] [5] [6] [7] , equation of 
motion for a domain wall without overhangs can be written in a following way 



where Z{x.,t) denotes the interface position; 7 and F are the domain wall 
mobility and stiffness, respectively; h is the external driving force. For a fer- 
romagnet h = jisHM , where H is the external magnetic flcld and Af is the 
magnetization. Finally rj denotes the random force generated by the impurities. 
For broad domain walls F ^ J/{a^~^l), where D denotes the dimensionality 
of the wall. For narrow walls F depends in general on J, T, and the disorder 
strength in a complicated way [8]. 

The random forces rj{x, Z) generated by pinning centers obey the Gaussian 
statistics and have short range correlations: 



Here 5i(x) denotes a delta- function smeared out over a distance /. The initial 
correlator Aq{Z)) is an even analytical function of Z, which decays to zero 



IdZ 
^~dt 



FV^Z + h + rj{x, Z) 



(1) 




(2) 
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over a finite distance I and has maximum ai z = 0. 

In the following we assume that the disorder is weak, i.e. that F ^ rjl. Under 
this assumption the interface is essentially flat on length scales L <^ Lp (see 

[6]), where Lp^ll-ij ^ / is the so-called Larkin length. On larger 



scales the wall adapts to the disorder and gets pinned for driving fields h ~ 

D/(4-D) 



hp with hp ^ TIL^ = rji^j ^ 77 for the pinning threshold. The 

transverse displacement w{L) of the DW at the scale L is determined by 
the roughness exponent ^(-^) oc ((^ < 1). The characteristic time of 
an avalanche t{L) at the scale L is determined by the dynamic exponent z: 
t{L) oc L\ 

If h exceeds hp, the wall starts to move. Close to the dcpinning transition 
the velocity, which can be considered as an order parameter of the transition, 
vanishes according to a power law 

, h> hp. (3) 




Here we introduced the characteristic velocity scale Vp = 'jhp. The moving 
domain wall conserves the rough (fractal) structure at length scales between 
Lp and a new scale established by small velocity ^„ = Lp{vp/v)~ , where the 
pinning force dominates. At larger scales the bumps on the DW are healed 
and it becomes smooth. 

The functional renormalization group calculations [6] [7] [9] show that 1 < z < 
2, i.e. the dynamics close to the depinning transition is super-diffusive, reflect- 
ing the rapid motion of the object after the maximum of the potential has 
been overcome. The critical exponents satisfy new scaling relations [6] 

.= ^ = ^>^, (4) 



These exponents were calculated first to order e = {A — D) in [6] and recently 
to order [9]: C = f (1 + 0.14331e) and z = 2 - | - 0.04321e2 [9]. 
In the opposite regime L ^ ^ the problem is essentially linear and Z can be 
replaced by vt in the argument of rj{'x, Z). This can be seen qualitatively as 
follows: On the time scale t the domain wall advances on average by an amount 
vt. Randomly distributed pinning centers will lead to a local distortion which, 
according to dynamical scaling, spreads over a region L{t) ^ Lp {tuptf^". The 
local retardation or advancement of the object due to the fluctuation in the 
density of the pinning centers scales as Z{t) ^ l{L{t) / LpY liupt)^/^, where 
uOp = Vp/l is the basic frequency at the Larkin length. Since C < -2, Z{t) grows 
more slowly than vt. Thus on time scales t > t^ = uj~^ {vp/vY^^^^^^ and length 
scales L > = L{ty) the non-linearities in the argument of ri{'x,vt + Z) can 
be neglected and the linearized theory applies. In this case ^(x, t) is replaced 
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by vt in the argument of the random forces. Random forces act then as a 
thermal noise with the temperature ~ v^^. 

After integrating out the interface fluctuations on the length scales L ^ ^y, 
the effective equation of motion for the interface profile Z(x, t) = (^(x, t))^^ 
where ()^^ denotes the spatial and time average over scales ^„ and ty, respec- 
tively, on large scales coincides formally with (1) with some values substituted 

by their renormalized counterparts: 7 — >• 7eff = 7 (^j^^ and h ^ h — hp. 
The influence of the random field on these length scales is negligible . 



2.2 Thermally activated creep of domain walls 



Next we consider the influence of thermal noise on the interface motion. Ther- 
mal fluctuations lead to a no-zero velocity of the domain wall even at small 
forces. All changes in h are still assumed to be made adiabatically. The right- 
hand side of equation of motion (1) must be complemented by thermal noise 
({x, Z, t) normalized as usual: 

{r){x, Z, t)r]{x', Z', t'))T = 2-5(x - x')d{Z - Z')d{t - t') . (5) 



In the absence of external driving force the typical free energy fluctuations of 
the length scale L are of the order F{L) = Tp{L/Lp)^ where x = -D — 2 -|- 2(^. 
The energy barriers between different metastable states scale have the same 
order of magnitude. When the driving force h switches on, it changes barriers 
between neighboring metastable states at the scale L by the value —hw{L) = 
—hLp{L/LpY. Thus, the total energy barrier is: 

EBiL,h) ^ F{L) - hL'^wniL) ^Tp^{l - \ (6) 



where we introduced the force length scale — Lp (^^^ A schematic 

graph of Eb{L, h) vs. h is shown in Figure 1, it has a maximum at L — Lh 
and vanishes for L — L^. The maximal height of the barrier is: 



It must be overcome by thermally activated hopping to initiate motion. The 
creep velocity of the domain wall follows from Vcreep ~ w {Lh) /r^Lh). According 
to the Arrhenius law, the hopping time is r ~ a;~^e^^('*)/^. Thus, we obtain 
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Fig. 1. Left: Energy barrier as a function of the length scale I for a given driving 
force density h. Right: The lines hoj and hr as explained in the text. 



Vcreep{h) ~ ^(i) ^ exp[— (^j ] . We have omitted a prcfactor, which is 



beyond our accuracy. This formula is valid for T ^ Esih) and was found first 
by loffe and Vinokur[15] (see also [16]). In the opposite case T 3> Esih) we 
expect a linear relation between the driving force and the velocity: v ~ 7/1. 
The border line between the two cases i.e. the inflection point of the curve 
v{h), T EB{h), defines a temperature dependent force hx (compare Figure 



Note that the creep formula is valid only for h ^ hx- 

Let us now consider the influence of thermal fluctuations on the depinning 
transition. At h < hp and T — the velocity is zero, but one has to expect 
that as soon as thermal fluctuations are switched on, the velocity will become 
flnite. Scaling theory predicts in this case an Ansatz [18,19] (generalizing (3)) 



with $(a;) const, for x and $(a;) ~ x^^^ for a; ^ 1, such that v{hp, T) ~ 

is a new exponent which still has to be determined . 
It is worthwhile to note that relevant thermal fluctuations, which unpin the 
domain wall act on scales of the order of Lp -C ^ as was flrst indicated by A. 
Middleton [19]. At the critical point h = hp essentially only barriers on the 
scale L Lp are left as we saw earlier. It is therefore sufficient to consider 
only this length scale. A detailed analysis of the form of the effective potential 
on this scale gives r = 3/2. 




1-) 




(8) 




(9) 
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3 Adiabatic motion of a single domain wall driven by ac force 

3. 1 Dynamics of a rectilinear domain wall at zero temperature 

In this section we describe the motion of a DW, rectihnear at large scale [20]. 
The reason of the motion is the alternating driving force, which we assume 
to have a simple harmonic shape: h{t) — ho sin cut. In previous section we 
demonstrated that the DW roughness can be ignored on a time scale t > t^ 
and the length scale L > L^. Thus, the bending of the DW can be neglected if 
Loty <^ 1. Simultaneously this condition means that the process is adiabatic, 
i.e. that the velocity at any moment of time with high precision is equal to 
its stationary value corresponding to the value of the driving force h{t) at the 
same moment of time. The rectilinear (or plane) DW can be characterized by 
one coordinate only. We denote it Z{t). In adiabatic approximation it satisfies 
an obvious equation: 



where the scaling function f{x) behaves asymptotically as at small x and 
as X at large x. Instead of integrating it over time, we integrate it over field 
by the following change of coordinate: dt = — r^= Thus, we find a following 

ul^^/h^-h'^ 

expression for Z as function of the field h: 



This equation is correct for h > hp] smaller positive values of h the DW 
does not move and Z = Zq = const. Therefore, the magnetization remains 
constant until the amplitude of the driving field reaches the value hp. This 
value marks the first dynamical phase transition: the appearance of the hys- 
teresis loop. Now starting from Zq = —L at /i = 0, where L is the half of size 
of a rectangular sample, and applying positive magnetic field, one can increase 
Z until it takes its maximum value 



This equation is valid as long as the right-hand side of equation (12) is less 
than L. When it becomes larger, the maximal value of Z obviously remains 
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(11) 




(12) 
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Fig. 2. Shapes of the hysteresis loops and dynamical phase transitions 

equal to L. The value ht2 of the amplitude h^, at which Z^ax becomes equal 
to L determines the second dynamical phase transition: the complete rever- 
sal of magnetization. The hysteresis loop becomes symmetric with respect to 
inversion: h —h, M —M [Z — >■ — Z) (see Figure 2b). At smaller val- 
ues of Kq between and hi\ the reversal of the magnetization is incomplete, 
the hysteresis loop is symmetric with respect to reflection h — > —h, but not 
to inversion of magnetization (see Figure 2a). Finally, when Z reaches the 
value L for a quarter of period or less, i.e. at h < Hq, in a range of fields 
hr < h < ho the magnetization reaches the saturation (in a little more realis- 
tic model it becomes a single- valued function of the field) . We call these parts 
of the magnetization graph whiskers (see Figure 2c). The value hta of ho, at 
which the whiskers first appear determines the third dynamical transition. It 
can be found from the obvious equation: Z{ht3, hts, —L) = L, where we have 
used notations of equation (11) . Note that similar equation is valid for ht2- 
Z{ht2,ht2,—L) — L/2 From these equation we find a scahng relation for the 
phase transitions fields: 



It is accepted to characterize the hysteresis loop by the coercive force h^ i.e. 
by a value of the driving force at which the magnetization vanishes (Z = 0): 
Z{hci ho, —L) = 0. Earlier we have introduced the saturation field hr ■ Equa- 
tion for it is: Z{hr,ho,—L) — L. Equation (11) implies that these two fields 
divided by hp obey the scaling equations depending on 2 dimensionless ar- 
guments u = ^ and V — The reader can find details in the original 
article [20]. We should mention the asymptotic behavior at very big u or v: 
ht2^ ^ ; ht3 ^2ht2; /ic ~ \Jht2 {2ho - ht2); K ^ ^ht^ {2ho - hn). The criti- 
cal asymptotics near the mobility threshold (f — 1 <C 1) is more complicated. 
The details can be found in the same article [20]. 

Numerical MC simulation of the 2-dimensional Ising model with random bonds 




(13) 



where F{x) is a function inverse to G{y) — f{y' — 1) 
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and the Glauber dynamics [20] demonstrated that the DWs, initially numer- 
ous, quickly merge to a few ones with linear size of the same order as the 
size of the sample. This stage of the magnetization reversal is the longest and 
determines dynamics in total. Thus, our results are qualitatively correct for 
larger systems at the last and longest stage of the hysteresis process. The same 
numerical simulation has reproduced all types of the hysteresis loop predicted 
by theory of the single DW hysteresis. However, the quantitative details may 
be different. Moreover, since theory of multidomain hysteresis does not yet 
exist, there is no reliable estimate of the time at which only few domains 
remain. 



3.2 Motion of a domain wall at finite temperature 



As it was shown in subsection 2.2, the principal difference of the dynamics at 
finite temperature is that the DW can start to move at an arbitrarily weak 
driving force due to thermal activation processes. Thus, strictly speaking, 
there is no more mobility threshold for the constant driving force. At low 
temperature the threshold will be smeared. Besides the finite tcmpcratTire 
establishes new scales of the length Lt, time tr oc L^, force Ht and activation 
energy EB{T,h). In adiabatic regime it leads to appearance of an effective, 
temperature dependent threshold for the ac driving force. In this section we 
follow the work [21] 

If the driving force is alternating in time with frequency uj and amplitude 
ho, the barriers at the scale L, for which ujT{L,h) > 1, where T{L,h) = 
Toexp Eb{L, h) /T is the relaxation time for such a fluctuation, can not be 
overcome during one cycle of oscillation. Thus, the global motion of the DW 
may be initiated only when the maximum of value u!t{L, h) over L becomes 
of the order of one. From the condition uJTraax{h) = 1 we find a new, frequency 
and temperature dependent driving field /i^^, which plays the role of dynamic 
threshold. It obeys a following equation: 



hp 




(14) 



where A = — ln(a;ro) (reminder: tq = lJq^ is a microscopic hopping time: we 
assume cutq <S 1). At h < h^, there is no macroscopic motion of the wall;, 
though still its segments transfer between different mctastable states with 
avalanches whose development time is less than 27r/w and the length scale 

is less than Lomega = Lp In ^) ■■ This process gives rise to dissipation. 
Drift of the wall as a whole starts at h^ > h^. At uotq ^ 1, /i^ < hx- Note 
that equation (14) determines /i^^ as a monotonously decreasing function of 
temperature accepting the value /ip at T = 0. At a fixed temperature h^ is 
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a monotonously decreasing function of frequency. Thus, DW subject to the 
ac driving force either remains in rest at < h^, or moves due to thermal 
activation (creep regime) at h^^ < hg < hr, or it moves in the sliding regime 
at /i > /iT- At low temperatures these three regimes overlap with the critical 
behavior near h — hp. A schematic phase diagram of the DW motion in 
variables T — h is shown on the r.h.s of Figure 1. 

Having derived the velocity of DW motion as a function of the driving field 
at fixed temperature and frequency, we can follow the hysteresis curve using 
the adiabatic equation of motion similar to that at zero temperature (10), but 
containing new parameters T and hi^: 



where f{x,y) is a dimensionless function of dimensionless arguments which is 
equal zero at a; < and 1 at a; or y infinite. 

The shape of hysteresis loops at finite temperature is rather similar to their 
shape at zero temperature. In particular, all 3 dynamic phase transitions de- 
scribed in subsection 3.1 proceed at finite temperature as well, but the role of 
the threshold driving force plays the field h^^. At the first transition the ampli- 
tude Hq coincides with h^^. At the second and the third transition amplitudes 
the same values are determined by equations: 



Equations (??,15) display the scaling similar to that at zero temperature: the 
dimensionless ratios htn/ h^, are defined by the same function of dimensionless 
variable {n — l)ujL/{2'yhJ). 



4 Non-adiabatic motion of a wall driven by an ac field 

J,..! Zero temperature 

We again consider the motion caused by an ac driving force h{t) = ho sin {cut) 
in eq. (1) with a frequency uj <^ cUp — "yhp/l, but we will not assume adiabatic 
condition. We start with zero temperature. 

The finite frequency lv of the driving force acts as an infrared cutoff for the 
propagation of perturbations generated by pinning centers. As it follows from 




(15) 




(16) 
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Fig. 3. Left: Schematic frequency-field diagram for the depinning in an ac external 
field (with ho > hp): For < uj <^ u)p the depinning transition is smeared but traces 
of the cj = transition are seen in the frequency dependency of the velocity at 
h = hp. This feature disappears for a; » Wp. Right: Velocity hysteresis of a D = 1 
dimensional interface in a random environment. 

(1) these perturbations can propagate during one cycle of the external force 
up to the (renormalized) diffusion length = Lp{'jT /toLpY^^ = Lp{ujp/ uj)^/^ . 
If < Lp, (i.e. uj > Up) there is no renormalization and z has to be replaced 
by 2. During one cycle of the ac-drive, perturbations generated by local pin- 
ning centers affect the configuration of the domain wall only up to a scale L^^. 
If this scale is less than Lp, the resulting curvature force TIL~'^ is always larger 
than the pinning force and there is no pinning anymore. 
In the opposite case > Lp (i.e. u < Up), the pinning forces compensate the 
curvature forces at length scales larger than Lp. As a result of the adaption of 
the domain wall to the disorder, pinning forces are renormalized. This renor- 
malization is truncated at L^^. Contrary to the adiabatic limit — » 0, there is 
no sharp depinning transition at > 0. Indeed, a necessary condition for the 
existence of a sharp transition in the adiabatic case was the requirement that 
the fluctuations of the depinning threshold in a correlated volume of linear 
size ^, 6hp ^ /ip(Lp/^„)(^+'^)/^, are smaller than {h — hp), i.e., {D + ()u > 2 
[17]. For uj > the correlated volume has a maximal size and hence the 
fluctuations 6hp are given by 

6hp ^ hp {Lp/Lj^^^^/' = hp (c./c.p)(^+^)/(^^) . (17) 



Thus, different parts of the domain wall have different depinning thresh- 
olds - the depinning transition is smeared in the interval 6hp. For a better 
understanding of the velocity hysteresis we consider coupling between the 
different segments of the domain wall with the average lateral size L^. Ap- 
proaching the depinning transition from sufficiently large fields, h^ ^ hp (and 
UJ <^ Up), one observes the critical behavior of the adiabatic C3iSG clS long as 
^ L^. The equality C,^ ~ L^^ defines a field hco signaling a crossover to 
an inner critical region where singularities are truncated by L^^. Note that 
hco — hp = hp{u/up)^^^''^'^ > 6hp (cf. Figure 3). 

A new physical phenomenon occurring in the non-adiabatic regime is the hys- 
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teresis of velocity. It is shown on the r.h.s. of Figure 3, which was obtained by 
numerical simulation of one- dimensional version of equation (1). The physical 
reason for this effect is the hampering of the DW by the pinning centers. It 
stops the average motion when the force becomes small, but not zero. Then the 
elastic forces tend to return too stretched intervals of the DW and submit it 
an average velocity opposite to the applied force, as is seen on the same figure. 
This consideration explains the clockwise circulation on the central hysteresis 
loop and appearance of secondary loops. The reader is referred for details can 
to reference [24]. 



5 Experiments 

In this section we review several experimental works relevant to our topic. We 
start with the work by Budde et al. [25]. They studied the first order phase 
transition between a two-dimensional gas and two-dimensional solid in the first 
adsorbed monolayer of a noble gas (Ar, Kr, Xe) on the face (100) of NaCl. 
Adsorbed layer was in equilibrium with the 3d gas of the same atoms. Varying 
temperature linearly with time back and forth, they observed the hysteresis 
of the adsorbate density. The deviation of the temperature from transition 
point plays in this experiment the role of the driving force. They have found 
that the width of the hysteresis loop, which corresponds to the coercive force 
he in magnetic system, scales as {NTY^^ , where AT is the amplitude of the 
temperature oscillations. This result agrees with theoretical predictions made 
in the work [20] (see the end of subsection 3.1). This result suggests that the 
experimental parameters corresponded to a regime of the force much larger 
than the threshold value, so that the pinning force was negligible. 
In the works by Kleemann and coworkers [26,27] the authors studied the di- 
electric spectra in ferroelectric single crystals Sro.6i-xBao.39Nb206Cex^ in the 
vicinity of its transition temperature. They have found different behavior of 
the electric susceptibility x vs. frequency. At low frequencies they observed 
a power-like behavior x oc a;'"' with 7 between 0.2-0.7. Such a behavior the 
authors ascribe to the creep motion with a broad distribution of the acti- 
vation energy. At higher frequencies they observed the logarithmic behavior 
of X, which they treat as the reversible relaxation of the DW segments. The 
crossover between these two regimes they attribute to the dynamical phase 
transition aX, h — as it is predicted in the work [21] (see subsection 3.2). 
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